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ABSTRACT 



We present theoretical evidence that the exclusive/inclusive ratio of semileptonic A b -decays exceeds 
that of semileptonic B-decays where the experimental exclusive/inclusive ratio amounts to about 
^ 66%. We start from the observation that the spectator quark model provides a lower bound on the 
^ ■ leading order Isgur-Wise function of the A& — ► A c transition in terms of the corresponding B — ► 
D, D* mesonic Isgur-Wise function. Using experimental data for the B — » D, D* mesonic Isgur- 
Wise functions this bound is established. Applying a Bethe-Salpeter model including spectator 
quark interactions and a QCD sum rule estimate of the A& — > A c transition form factor which 
satisfy the spectator quark model bound we predict the exclusive/inclusive ratio of semileptonic 
At decay rates to lie in a range between 0.81 and 0.92. We also provide an upper bound on the 
baryonic Isgur-Wise function which is determined from the requirement that the exclusive rate 
should not exceed the inclusive rate. 
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I. INTRODUCTION 



In mesonic semileptonic b — > c transitions, the exclusive transitions to the ground state S'-wave 
mesons B —>■ D,D* make up approximately 66% of the total semileptonic B — > X c rate It 
would then be interesting to know what the corresponding semileptonic rate ratio T^ b ^,^ c /T^ b ^x c 
(termed Re in the following) is in semileptonic A^-decays. This is an important experimental 
issue since a knowledge of this ratio would greatly facilitate the analysis of semileptonic A^-decays. 
For example, if the semileptonic A^-decays were dominated by the quasi-elastic exclusive channel 
Ab — > A c + l~ + Vi this would be of considerable help in the kinematical reconstruction of 
their decays in as much as the A c baryon is easy to detect via its decay mode A c — > pK~7i + . 
Unfortunately nothing is known experimentally about this ratio yet. 

In this paper we attempt to address the problem of a determining the exclusive/inclusive 
ratio Re in semileptonic A^-decays from a theoretical point of view by consulting some model 
calculations which we critically scrutinize. We also attempt to extrapolate from the experimentally 
known results in the meson sector to the baryon sector. 

As concerns the inclusive semileptonic rates of bottom mesons and bottom baryons one is 
now reasonably confident that they can be reliably calculated using the usual operator product 
expansion within HQET. The leading term in the OPE is given by the free heavy quark decay 
rate which clearly is the same for baryons and mesons. Radiative corrections to the free quark 
decay rate are quite large but again are identical for mesons and baryons. Differences in the 
inclusive semileptonic rates of mesons and the Ab baryon set in only at 0(1/ ml). They affect the 
mesonic and A& rates differently since there is no chromomagnetic Oil/ml) correction in the A& 
case. However, since the chromomagnetic term contributes only at the 3.7% level, the difference 
in the inclusive semileptonic rates for mesons and baryons is predicted to be quite small. 

A much more difficult task is to get a reliable theoretical handle on the quasi-elastic exclusive 
semileptonic Ab — > A c rate. There exist a number of theoretical calculations on the exclusive 
decay A& — > A c + l~ + V\ using various model assumptions. They are of no great help since their 
predicted rate values may differ by factors of up to three and it is not easy to judge the reliability 
of the various model assumptions that enter the calculation. Ideally one would like to have model 
calculations that are valid both in the heavy meson and the heavy baryon sector. If these model 
calculations give sensible results in the heavy meson sector, where they can be checked against 
data, one would have more confidence in their predictions for the heavy baryon sector. 

The paper is structured as follows. In Sec. II we take a first look at the leading order rate formula 
for the exclusive semileptonic decays of B mesons and A& baryons in order to get a semiquantitative 
handle on the relative size of their rates. The analysis is refined in Sec. Ill for A& baryons where 
l/m& and corrections and renormalization effects are included. In Sec. IV we recapitulate 

the calculation of the semileptonic inclusive decay rates. The results of the Sections III and IV 
are brought together in Sec.V where we discuss the exclusive/inclusive ratio Re of semileptonic 
Ab decays. We present numerical results on the exclusive/inclusive ratio for various models and 
give our best estimate of this ratio. In Sec. VI we classify the possible non-exclusive final states 
that will have to fill the gap between the exclusive and inclusive rates in semileptonic A& decays. 
Sec. VII, finally, contains our conclusions. 

II. THE HEAVY QUARK LIMIT 

For a quick first appraisal of the question of how the exclusive semileptonic decays of mesons 
and baryons are related we turn to the heavy quark limit and list the leading order semileptonic 
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rate formulae for the B — > (D + D*) and A5 — > A c transitions. One has |]|| 

^^"^ = ^%^r 3 v^T(3^(l + r 2 ) - 2r(2^ 2 + 1)) { F ™ ^ 1 * ) , (2.1) 

where r = M 2 /M 1 and = (M 2 + Mf — q 2 )/2MiM 2 . F meson (u) and -Fbaryon(w) are the leading 
order Isgur-Wise transition form factors for the B — > D,D* and — > A c transitions, respectively. 
Throughout the paper we refer to Mi and M 2 as the masses of the initial and final particles in the 
semileptonic decay process. 

The free heavy quark decay rate (or leading order parton model rate) which we need later on 
is simply obtained by replacing the particle masses in ( |2.1| ) by the corresponding quark masses 
and setting the curly bracket in fl2.1|) to one, i.e. by taking the current coupling in the A& — > A c 
case to be point-like. We shall encounter the integrated parton model rate for the A& — > A c case 
again in Sec. IV. Finally, in the heavy quark limit one has to determine the final meson mass M 2 
by taking the weighted average M D = 1/4(M D + 3M D *) = 1.973 GeV with M D = 1.869 GeV 
and M D * = 2.010 GeV. For the pseudoscalar bottom mass we take M B = 5.279 GeV. For the 
Ag-baryon masses we use M^ b = 5.624 GeV and Ma c = 2.285 GeV. 



When trying to compare the two rates in ( |2.1| ) one identifies two main determining factors 



which counteract each other. On the one hand one has the form factor expressions in the curly 
bracket which tend to enhance the mesonic rate due to the multiplicative factor {u + l)/2 in the 
mesonic case. Also, according to common prejudice the baryon form factor falls off more rapidly 
than the mesonic form factor. On the other hand one has the overall mass factor Mfr 3 which 
enhances the baryonic rate because Mfr 3 = 214.16 GeV 5 and Mfr 3 = 377.37 GeV 5 in the mesonic 
and baryonic cases, respectively. 

It is evident that the choice of mesonic and baryonic Isgur-Wise functions plays a crucial role 
when comparing the two rates. As has been emphasized before, there exist some experimental 
knowledge on the mesonic Isgur-Wise function but nothing is known experimentally about the 
baryonic Isgur-Wise function yet |. 

For quick reference it is sometimes convenient to characterize the fall-off behaviour of the 
Isgur-Wise functions by expanding it around the zero recoil point where one has the zero recoil 
normalization condition -F(l) = 1. Keeping terms up to second order in this expansion one has 

F(u) = F(l)[l - P \uj - 1) + c (w - l) 2 + ... ], (2.2) 

where the coefficients p 2 and c are called the slope parameter and the convexity parameter, re- 
spectively. The slope parameter is frequently used to characterize the fall-off behaviour of the 
Isgur-Wise function. The expansion ( |2~2"1) is useful if one studies the physics close to zero threshold 
but may give misleading results when calculating rates because the spectral weight function mul- 
tiplying the form factor functions is essentially determined by the square root factor sJuS 1 — 1 in 
( |2.1| ) and is therefore strongly weighted towards the end of the spectrum. It goes without saying 
that the slope and convexity parameters are in general different for the mesonic and baryonic 
Isgur-Wise functions. 



3 The only available experimental result is from a preprint version of a DELPHI-analysis Q. This paper 
quotes a value of p 2 = 1.81^Q'g^ ± 0.32 for the slope of the baryonic Isgur-Wise function. However, since 
this paper has never been published, we shall not use their result in our analysis. 
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FIG. 1. Leading order baryon form factor -Fbaryon^) m the dynamic Bethe-Salpeter model of Q 
(denoted by C( w ) m this figure), a) Noninteracting light quarks (long-dashed line) and b) interacting 
light quarks with the range parameter As = 500 MeV (solid line) and Ab = 355 MeV (short-dashed line). 



In order to proceed with our first appraisal of the magnitude of the exclusive mesonic and 
baryonic semileptonic rates we appeal to the spectator quark model where the mesonic and baryonic 
Isgur-Wise functions become related to one another |2|,|4j]. In the spectator quark model one has 

-fbaryon(^) = ^ | -^meson (^) | • (2-3) 

Explicit calculations show that the baryonic form factor is considerably underestimated by the 
spectator relation ( ^3|) . Nevertheless, the spectator relation ( |2.3|) may still serve as an effective 
lower bound on the baryonic form factor. 

The physical picture behind the spectator quark model relation is quite simple. In the heavy 
baryon case there are two light spectator quarks that need to be accelerated in the current transition 
compared to the one spectator quark in the heavy meson transition. Thus the baryonic form factor 
is determined in terms of the square of the mesonic form factor. The factor is a relativistic 

factor which insures the correct threshold behaviour of the baryonic form factor in the crossed 
e + e~-channel 

In [[|] the relation between heavy meson and heavy baryon form factors was investigated in 
the context of a dynamical Bethe-Salpeter (BS) model. The above spectator quark model relation 
( p.3|) in fact emerges when the interaction between the light quarks in the heavy baryon is switched 
off in the BS-interaction kernel. In the more realistic situation when the light quarks interact with 
each other, the heavy baryon form factor becomes flatter, i.e. the spectator quark model form 
factor may be used to bound the heavy baryon form factor from below. In Fig.l we reproduce 
from 0] the u dependence of the spectator quark model form factor and that of two representative 
form factors with the interaction between the light quarks included. The starting point in is a 
mesonic form factor with a slope of p^eson = 1 which, according to (|2.3|) , leads to a spectator form 
factor with a slope of Pb aryon = 1-5- The spectator form factor is the lowest form factor shown in 
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Fig.l. The interaction between the light quarks was introduced through a harmonic oscillator type 
kernel in the BS-equation. The two upper form factor curves in Fig.l correspond to two different 
choices of the oscillator strength with which the light quarks interact or, equivalently, correspond 
to two different choices of the size parameter in the oscillator wave function. The interaction type 
form factors in Fig.l have slopes of Pbaryon = 0-81 (solid line) and 0.97 (short-dashed line) Q. 
They are considerably flatter than the spectator quark model form factor. 

We shall now calculate exclusive rates for mesonic and baryonic transitions according to ( |2.1| ) 
using the spectator model relation ( |2.3| ). According to what was said before, the baryonic rate 
calculated in this way must subsequently be adjusted upward according to the analysis of [f|]. We 
went to considerable lengths in explaining the results of [|J because we want to emphasize that 
the outcome of the baryonic rate estimate using the spectator quark model relation (|2.3|) must be 
viewed as providing only lower bounds on the true quasi-elastic baryonic rate. For the mesonic form 
factor we use the world average of the slope p^eson obtained by combining results from B — > D* 
and B — > D transitions, p^eson = 0-^0 0- Using V c b = 0.038, a linear meson form factor with the 
above slope, a baryonic form factor according to the spectator relation ( |2.3j ) and the rate formulae 
Q2.1| ) one obtains r meson = 5.30-10 10 s _1 and Tbaryon = 5.04-10 10 s~ 1 . As has been emphasized before 
the baryonic rate Tbaryon = 5.04 • lO 10 ^" 1 has to be adjusted upward in the more realistic situation 
of interacting light quarks. Looking at the model calculation M for guidance, the increment in rate 
going from noninteracting (spectator) to interacting quarks is 1.28 and 1.37, respectively, for the 
two choices of oscillator strengths analyzed in ||. Adjusting the above baryonic rate accordingly 
our leading order estimate of the baryonic rate is thus Tbaryon = (6.45 -j- 6.90) • 10 10 s _1 . Starting 
from a mesonic exclusive/inclusive ratio of 66% and assuming equal inclusive semileptonic rates 
for bottom baryons and mesons, which is sufficiently accurate for our semiquantitative calculation, 
our estimate for the exclusive/inclusive ratio in semileptonic A b decays is Re = (80 — 86)%. This 
is considerably larger than the mesonic exclusive/inclusive ratio Re ~ 66%. 

Up to this point our semiquantitive analysis was done to leading order in HQET. How would 
finite mass effects affect our previous conclusions? One way of improving the previous analysis in 
the meson sector is to insert physical masses in the rate expression ( |2.1| ) , thereby including part of 
the 1/mQ-corrections to (|2.1|) . To do this we need to disentangle the B — > D and B — > D* rates 
in (Q. One has §: 

dT ^ D) = ^^^' ^(l + r) V - lplFmesonMI 2 (2.4) 

and 

dT ( V ^ = ^^r 3 V^l(u + 1) [(1 - r) 2 ( W + 1) + 4u,(l - 2ur + r 2 )j |F mcson (^)| 2 . 

(2.5) 

Using now the world average of (p 2 I1 eson)' B_+D = 0-66 and (Pmeson) B ^ D * = 0-71 with linear form 
factors and taking physical D and D* masses one finds Tb->d — 1-39 • 10 10 s^ 1 and r^^^. = 
3.90 • 10 10 s _1 giving a total mesonic rate of T B ^ D+D ,* = 5.30 • 10 10 s~ 1 . It fully agrees with the 
above result and therefore leaves the aforegoing conclusions intact. 



4 The CLEO Coll. also attempted a linear plus quadratic fit to the data, but the data was not good 
enough to determine the convexity parameter c of the meson form factor with any accuracy. 
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Continuing with our discussion on the contributions of nonleading effects in the 1/mQ-expansion 
we now turn to results of some model calculations in order to find out how nonleading effects may 
affect the above conclusions. In the mesonic sector Neubert and Rieckert || analyzed an infinite 
momentum frame model and found that 0(l/mg) effects raise the B — > D and B — > D* rates by 
15.7% and 0.5%, respectively, resulting in a rise of 4.4% for the total D + D* rate. Using a similar 
infinite momentum frame model Konig et al. find that the 0(l/rriQ) effects raise the semileptonic 
Ab — > A c rate by 3% j?J which is quite close to the 4.4% found in [BJ in the bottom meson case. 
Judging from these model calculations our leading order comparison of the mesonic and baryonic 
rates and the conclusions drawn from it do not seem to be much affected by 0(1/ 'w,q) corrections. 

There also exist estimates of 0(l/mg) corrections in the literature. Faustov and Galkin et 
al. use a relativistic quark model based on the quasipotential approach ||. They quote exclu- 
sive/inclusive branching ratios of (13.5 + 3.3 — 1.4)% and (39.1 + 6.5 — 3.9)% for semileptonic 
B — > D and B — > D* rates, where the second and third numbers refer to the O{l/mo) and 
0(l/m 2 Q) corrections, respectively. The (9(l/mg) corrections in this model are considerably larger 
than in the infinite momentum frame models. Ivanov et al. investigated the role of finite mass 
effects in semileptonic A;, — > A c decays without taking recourse to the heavy mass expansion. They 
found an overall rate reduction of 9% relative to the infinite mass result ||. In the analysis of the 
present paper presented in Sec. Ill we obtain « +5% and ~ —7% for the 0(l/rriQ) and 0{l/mn) 
corrections to the A^ — > A c decays, respectively. From all these model calculations one learns 
that the 0{\jmQ) corrections tend to increase the rates whereas the 0(1/toq) corrections tend 
to decrease the rates, for both heavy meson and heavy baryon decays. Again, our leading order 
estimate of the relative size of the exclusive/inclusive ratios of bottom mesons and baryons is not 
likely to be affected much by including also (9(l/mg) effects. The same holds true for renormal- 
ization effects of the weak current which affect the bottom baryon and bottom meson amplitudes 
equally and therefore drop out in the ratio of exclusive semileptonic bottom meson and baryon 
decays. 

The conclusion drawn in this section on the predominance of the exclusive/inclusive ratio of 
semileptonic Af,-decays over that of semileptonic B-decays carries over to the more sophisticated 
analysis of the next sections where we include 1/mq and l/m| effects, and radiative corrections. 



III. EXCLUSIVE SEMILEPTONIC RATE A 6 -»• A c + l~ + V t 



It is most convenient to represent the differential decay rate in terms of the helicity amplitudes 
of the process. One has (we take leptons to be massless) [P,|I~0H 



ffiir+i^_i_ 1 i 2 +iffi r+i^_ioi 2 ). (3.i 



dY (A„ -> A c ) = G 2 F \V bc \ 2 g 2 Mfy/^T 
du 96vr 3 Mi 

The helicity amplitudes are in turn related to the invariant amplitudes of the process via 



'q 2 HY' A = pM 1 M 2 {u T 1) ((Mx ± M 2 )f? A ± M 2 (lo ± l)f 2 v ' A ± M x {u ± l)ft A 



Hi* = -2^M x M 2 (uj =p l)/r , (3.2) 

2 1 v 

where the invariant amplitudes are defined by 

< A c (v 2 )\jV\A b (vi) > = m c (w 2 )(/i / 7 m + f?v lfl + f£v2p)u b (vi), 

< A e (v 2 )\ J A \A b (vi) > = uMifiln + f 2 v lfl + fiv 2ll )^u b { Vl ). (3.3) 
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The total helicity amplitudes finally are given by 



A2AW 



(3.4) 



where the choice of of the relative minus sign between the vector and axial vector helicity am- 
plitudes reflects the (V — A) structure of the b — > c current transition. The are the 
helicity amplitudes for the vector (V) and axial vector (A) current induced transition in the decay 
1/2+ -> 1/2+ + W, 



of f— shell 

boson, respectively. 



with A2 and Xw being the helicities of the final state baryon and the W 

by parity. 



The remaining helicity amplitudes are related to the above helicity amplitudes 
One has 



H 



V,A 

— A2— Av 



±H 



V,A 
A2AH 



(3.5) 



It is well known that the complexity of the form factor structure exemplified by the set of six 



form factors /, 



V,A 



1, 2, 3) is considerably reduced in HQET. Working up to 0(1/ 'ttiq) in HQET 



and including also 0(a s ) corrections one finds [|| 

1 1 



A» = F(a;) + 



/ 2 » 



F(u) 



2Mi 2M 2 

I 1 X _ a s( m 



H 0*0 = F(u) (- 
tf{u>) = F(u) + 



f 2 A (uj) = F{u) 



M 2 uj + l 

1 1 
M 1 u + 1' 

( 1 1 



V (uj) + AF(lu] 



a Am) 



71 



vi(uj, X)F(u) 



-A 



7T 

a Am) 



71 



-V 3 [U 



V2Mi 2M 2 



/ \ —r- ^ / k UJ — 1 

r}(u) + AF(u) 



a s \m) 

7Y 



ai(u, X)F(uj) 



#{u>) = F(u) 



M 2 uj + 1 
Miu+l 



1 -r a s (m) , . 
A — — -a 2 (uj) 



71 

aJm) 



TV 



-a 3 {u 



(3.6) 



The 0(a s ) corrections to the form factors have been taken from 
0(a s ) vertex correction to the current-induced b 



11 1 . They result from the 



c transition 12 



The infrared singularity is 
regularized by the introduction of a fictitious gluon mass which is taken to be A = 0.2 GeV. At 
zero recoil, where the vertex correction is infrared finite, the renormalization is independent of 
the gluon mass regulator. However, away from zero recoil, the a s -correction functions vi(u,X) 
and ai(o>, A) depend on the gluon mass regulator. This introduces a certain amount of model 
dependence in the renormalization procedure. The above value of the gluon mass was chosen 
according to the expectation that the exchange of virtual gluons in the vertex correction should 
be cut-off at frequencies k° ~ 1/R with R ~ 1 fm being a typical hadronic scale. The argument of 
the a s coupling, m, is taken such that effects of higher order terms (a s \n(mf ) /m c )) n are minimized, 
m = 2m(,m c /(m( ) + m c ) ~ 2.31 GeV. In this way one avoids the use of the renormalization group 



np. 



improved summation of leading logarithms, which has been proven as inconsistent 

The binding energy of the A& is denoted by A which we take to be 0.6 GeV. The form fac- 
tor function 77(2;) results from the nonlocal contribution of the kinetic energy term of the 1/mq 
corrected HQET Lagrangian. It has been calculated in two different model approaches and has 
found to be negligibly small J?, 13]. Therefore, we can safely drop its contribution in the follow- 
ing. By neglecting the form factor r](x) in the O(l/mo) result (|3.6| ) the differential rate ( |3.1|) 
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is proportional to the square of the leading order Isgur-Wise function F(u) with its zero recoil 
normalization F(l) = 1. In this way we can meaningfully compare our results with the leading 
order results of other model calculations as will be done in Sec.V. 

It is well known that 0(l/mq) corrections to the unit zero recoil normalization of Isgur-Wise 
functions can be substantial. For example, by evaluating zero recoil sum rules, the authors of Refs. 



14 obtain 



F(1) b -d = 0.98 ±0.07, 

F(1) B ^ D * = 0.91 ± 0.03. (3.7) 

in the mesonic case. In the A& — > A c case the zero recoil sum rule gives a bound on the zero recoil 
value of the sole remaining form factor function ff. The (unrenormalized) zero recoil sum rule 



leads to ff(l) < (1 - 0.165 /4/GeV 2 ) 1 / 2 |§. Using /i 2 = 0.5 GeV 2 [18,19], ff(l) must be smaller 



than 0.958. For defmiteness we take a value close to the upper bound 

A A (1) = 0.95. (3.8) 

This value is nicely corroborated by the finite heavy quark mass calculation of |J where one finds 
./Vd) 0.97. 

Nothing is known about the size of the 0(1/ ran) corrections to A& — > A c away from zero recoil, 
except that they can be parametrized in terms of ten new tu-dependent form factors and one new 
dimensionful constant [|16|], the magnitude and functional forms of which are not known. The lack 



of knowledge about the 0(l/mn) corrections away from zero recoil prevents us from their exact 
treatment. On the other hand the size of the 0(1/ nig) correction at zero recoil is a clear indication 
that the 0(1/ w/q) corrections cannot be neglected. We shall therefore adopt the following strategy. 
We smoothly extrapolate from the 0(1/ toq) information at zero recoil to the whole cu-range. The 
appropiate amplitude for this extrapolation is the axial vector current S-wave amplitude the zero 
recoil value of which is determined by the zero recoil sum rules. We thus multiply the axial vector 
current S'-wave amplitude everywhere by its zero recoil value fi(l) = 0.95. It is clear that the 
0(1/171%) corrections at zero recoil are exactly included in this approach. The lack of knowledge 
about the 0(1/ run) corrections to the other partial wave amplitudes leaves us no choice but to 
leave them untreated. With this in mind it is gratifying to note that the S'-wave contribution 
dominates the quasielastic rate. For example, using the standard form factor (|5.1|) with p\ = 0.75 
in a leading order calculation one finds that the S-wave contribution amounts to ~ 66% of the 
total semileptonic rate. 

In order to set up our procedure of how to incorporate the 0(1/ mq) corrections we define the 
relevant vector current (V) and axial vector current (A) partial wave amplitudes A^'$ in terms of 
the helicity amplitudes H x £\ w - L denotes the orbital angular momentum of the final state and 
S = ./current + Sa c is the sum of the final state spin angular momenta where J CU rrent = 1 in the zero 
lepton mass case that we are considering here. One has 



4i 

1 2 




A* =K\-Ht-J-Ht 

°2 V3 2° V3 2 1 





= \ hHt + \ hHt- (3-9) 
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By substituting invariant form factors according to ( |3.2|) one can verify the correct threshold 
behaviour of the partial wave amplitudes, i.e. AY!,AY 3 ~ (u — l) 1 ^ 2 , A^i ~ (u — 1)° and t4A ~ 

1 2 2 U 2 2 

According to the above strategy we now incorporate the OiA/mq) corrections by multiplying 
the S"-wave amplitude Afii by the C(l/mg) zero recoil correction ff(l) = 0.95. Thus we write 

<! - = <| + (tf (i) " (3-10) 

For the first term on the r.h.s. of ( |3.10| ) we substitute the 0(1/ uiq) result according to (|3.(j|). 
Contrary to this we use only the leading order result for A^i in the second term of fl3.10| ) since 

it is already being multiplied by the 0(l/m 2 n) factor (fi(l) — 1). Including also the A^ 3 partial 
wave amplitude the leading order expressions for the axial vector partial wave amplitudes read 

1 2 



u 2 



aJm] 



A± 



71 

1 2V2 

7T 



^M x M 2 (lo + l)F(cu) 
' (w - l)(M 2 a 2 (u;) + Mia 3 (a;)) 



( Ml - M 2 + 2^(1 + ^^ ai (w, A)) 



7T 



MiM 2 (a; + l)F(a; 
w- l)(M 2 a 2 (cj) + M 1 a 3 (cj)) 



( Ml _ M2 _ J^)(i + ^l ai {u, A)) 



7T 



(3.11] 



Putting everything together we arrive at the differential rate. One obtains 
dT (A b - A c ) G 2 F \V bc \ 2 q 2 M*V^T r v 2 y 2 A 2 A 2 



du 



3 q z 



( Ml - M 2 + 2^(1 + 2^pL 1 (u;, A)) 



■M 2 + 2^q 
aJrft 



71 



71 



co- l)(M 2 a 2 (u) +M 1 a 3 (u)) 



. (3.12) 



In Eq. Q3.12D the first line incorporates the 0(1) and 0(1/ tuq) contributions including the radiative 
corrections as specified by ( |3.2|) and by (|3.6| ). The second and third lines comprise the 0(l/mq) 
corrections (including radiative corrections) as described before. Since our aim is to compare our 
exclusive rate with the inclusive 0(a s ) rate written down in Sec. IV we have only retained radiative 
corrections up to 0(a s ) in the exclusive rate ( |3.12| ) for consistency reasons. 

For the sake of completeness we separately list the 0(1/ m^) zero recoil corrections for the 
longitudinal and transverse pieces of the axial vector contributions. They are needed for the 
transverse/longitudinal separation shown in Fig.2. One has 



A |2 

ol 



HU 2 + ((f?(l)f - l)i:(A^) 2 + X r-(ft (1) - • K 



Ht 2l ? - + ((A A (1)) 2 " l)^«i) 2 " ^(^(1) - ^ ■ A 



rAl 2 

2 



^2 

2 ' 



(3.13) 



When summing the two contributions ( |3.13| ) in the rate formula the A^i ■ A^3 interference con- 
tributions cancel out as is apparent in (|3.12|) . As explained before we shall use the leading order 
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FIG. 2. u spectrum of the exclusive decay rate and partial rates of longitudinal and transversal 



transitions calculated with the standard form factor ( |5.1|) using p B = 0.75. 



results ( |3.11|) for the second and third term in (|3.13|) since the factors ((/^(l)) 2 — 1) and (ff(l) — 1) 
multiplying them are already of OiXjm 2 ^). 

Our numerical evaluation of ( |3.12| ) is based on the standard form factor (|5.1| ) with p\ = 0.75 
using again Vj, c = 0.038. All parameters have been specified before. For the quasi-elastic rate we 
find r excl = 5.52 x 10 10 s _1 . The 0(l/m Q ) and 0{l/m 2 Q ) corrections amount to +5.2% and -6.6%. 
The renormalization of the heavy quark current decreases the exclusive rate by 8.8%. In Fig. 2 we 
show a plot of the tu-spectrum of the quasi-elastic rate where we separately show the transverse 
(Avk = ±1) and longitudinal contributions (\w = 0) including C(l/mg) and 0(a s ) corrections 
calculated according to (|3.13|) . The longitudinal rate dominates the spectrum except for a small 
region close to zero recoil. For the integrated rates we find Tf^ cl /T^ cl = 1.89. 



IV. INCLUSIVE SEMILEPTONIC RATE A 6 -► X c + l~ + F, 

To the leading order in the heavy mass expansion the inclusive rate is given by the free heavy 
quark decay rate which can be obtained from Eq. (|2.1|) using quark masses and setting the curly 
bracket equal to 1 . There are no 0(1/ tuq) corrections to this result. Mass corrections come 
in at the order 0(l/m^). In the case of A&-decay, where the light diquark system has spin 
0, the chromomagnetic contribution drops out and the mass corrections are determined by the 
nonperturbative kinetic energy parameter fi 2 alone. 

Including also the a s -correction in the free quark decay rate [17], one has (x = (m c /mb) 2 ) 



i- = r „ f i - l^ 9 (x)\ (i - A] , (4.i) 



7T J \ 2m 

where r is the lowest order (in a s ) free quark decay rate 
^ G 2 F \V bc \ 2 ml 



1927T 3 



-I (x) , I (x) = (1 -x 2 )(l -8x + x 2 ) - 12x 2 \nx, (4.2) 
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and the function g(x) is determined by the 0(a s ) radiative corrections including all mass correc- 
tions as calculated in |17| : 



g(x) = h(x)/I (x) , 
h(x) = -(1-x 2 ) ( 



25 239 25 9 \ , . , /' 1 ., 17 ., , ., ., 

— 3~ X + X ) + xin ( x ) ( 20 + 90x ~ % x + y x ) +x ln (a0(36 + ar* 



2 . /17 64 17 2 \ , ^ . 
+ (1 - x 2 ) ^— - — + — x 2 J ln(l - x) 

- (1 + 16x 2 + x 4 )[6Li 2 (x) - vr 2 ] 



4(1 + 30x 2 + x 4 ) ln(x) ln(l - s) 



32x 3/2 (l 



7T 



- 4Li 2 ( % /x) + 4Li 2 (-v^) - 2 ln(x) ln 



1- v^' 



(4.3) 



The 0(1/ m 2 ) corrections appear in the third factor of Eq.(OJ). For the value of the kinetic energy 
parameter // 2 we take [|T8|JT9f 



/i 2 = 0.5 GeV 2 , 



(4.4) 



where we assume equality of the kinetic energy parameter in the meson and baryon case. 

It is well known and evident from Eq. fl4.2f) that the inclusive decay rate depends rather strongly 
on the exact value of the 6-quark mass m& which is fraught with some uncertainties. We shall use 
the results of two recent theoretical analyses of the inclusive semileptonic decay rate. In |19| the 



value of mi, was determined from an analysis of T sum rules and 5-meson semileptonic widths: 



m 6 = 4.8GeV, m c — 1.325 GeV, 

where the charm quark mass was determined from the constraint 

1 1 



2m c 2m& 



M B - M D . 



(4.5) 



(4.6) 



The M b,d are the spin-averaged masses M B D = l/4(M BjD + 3M B * D *) as before. 



In [20] the inclusive semileptonic B decay rate was directly expressed in terms of the T(IS') 
meson mass instead of the b quark mass. The authors of P0[ obtained: 



lincl— Y 



fm r \ ! 
192tt 3 V 2 J 



0.533 



0.096e - 0.029e 2 - (0.28A 2 + 0.12Ai)/GeV 2 



(4.7) 



where e = 1 denotes the order of the expansion in mr- Mass corrections and radiative corrections 
are already taken into account. The parameters Ai and A 2 in Eq.(3~7) are connected with the more 
familiar /i 2 and jji 2 G parameters by /i 2 = — Ai and /Jq = 3A 2 = 0.36 GeV 2 . For the A& baryon we 
set A 2 = and assume equality of Ai in the meson and baryon before. 
For the semileptonic inclusive b — > c decay rate of the Af, we finally obtain r mcl 



6.50- lO 10 ,^ 1 



using the mass parameters from |19| and r incl ~ T = 6.23 • 10 10 s _1 using the evaluation of ||20|| . 
Again we have set He = 0.038. We mention that these two inclusive rate values include the 0(a s ) 
radiative corrections which lower the inclusive rates by about 11%. This value is not very far away 
from the 8.8% by which the exclusive rate gets lowered by the same radiative corrections. 
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V. THE EXCLUSIVE/INCLUSIVE A b -> A c RATIO 

In this section we determine the exclusive/inclusive ratio Re = r excl /T mcl in semileptonic 
Afc — > A c decays based on our estimates for the inclusive rates derived in Sec. IV and on various 
phenomenological models for the baryonic Isgur-Wise function F(u) entering in the exclusive 
differential rate Eq. ( |3.12| ). Of all the phenomenological models we shall mostly focus our attention 



on the sum rule calculation of 21 



We begin our discussion with the determination of the leading order A& — > A c Isgur-Wise 
function by the QCD sum rule method given in ||21|| . The shape of the Isgur-Wise function in [21| 



can be very well reproduced by an exponential representation of the form 

F(u) = exp (-(2p| - 1)^) , (5.1) 

which has the correct zero recoil normalization F(l) = 1 and a slope parameter given by p\. 
The convexity parameter in this representation (proportional to {u — l) 2 )) is given by cb = 
l/8(-l+4p| + 4p|) and is positive for p\ > 0.207 as in most model calculations. We refer to this 
representation of the Isgur-Wise function as the standard form. For the p\ parameter the authors 



of IflTJ find p B = 0.85 and p B = 0.65 using diagonal and nondiagonal sum rules, respectively. As 



an average of these two values one obtains 

p\ = 0.75. (5.2) 

Using the average value of p B , Vb c = 0.038, the standard representation of the Isgur-Wise function 
7i]) and the rate formula ( |3.12| ) from Sec. Ill one obtains the exclusive rate T excl = 5.52 • 10 10 s _1 . 

^incl _ c t(l . 1 nl0„-l 



From the inclusive rate calculated using the mass parameters given in |TJ r incl = 6.50 • 10 10 s 
one finds Re = 0.85 for the exclusive/inclusive ratio. Note that the He dependence drops out in 
this ratio. The values of the slope parameter p\ and the exclusive/inclusive ratio Re of the model 
of [21] as well as those of other phenomenological models have been collected together in Table I. 
Table I uses the larger value of the two inclusive reference rates discussed in Sec. IV based on the 
mass parameters of ||19|| . If one instead uses the inclusive rate of ^(J all Re- values in Table I have 



to be increased by 4.3%. Radiative corrections do not affect the exclusive/inclusive ratios listed in 
Table I very much since they lower both the exclusive and inclusive rates (see Sees. Ill and IV). 
If they were left out the exclusive/inclusive ratio Re would be reduced by w 2%. 

It is clear from Table I that the form factor calculated in the quark confinement model |23| is 



too flat to satisfy the bound r excl < r . All other models in Table I satisfy this upper bound. 
Translating the upper bound Re = 1 into a lower bound on the slope parameter p\ one obtains 

(p|)min = 0.36 (5.3) 

using the standard form factor function (|5.1|) and the inclusive rate calculated from the mass 



parameters in [jHJ 



An upper bound on the slope parameter ps can be obtained using the spectator model bound 
discussed in Sec. II, which reads 

P% < 2Pm ~ \ (5-4) 

The mesonic slope parameter p 2 M can be extracted from the exclusive semileptonic B decays 0. 
The values are 
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TABLE I. Predictions of the slope parameter p 2 B and the ratio Re = r excl /r mcl calculated in different 
models. 

MODEL ISGUR-WISE FUNCTION p% r cxcl /T incl 



quark confinement model §1 in(a;+V^ 2 ^_L) 0.33 1.02 

Vu — 1 



QCD sum rules 13 J ^ exp (-0.8^ J 0.65 0.89 

QCD sum rules (p 2 = 0.75) ^ exp (-{2p 2 - 1)^) 0.75 0.85 

Simple quark model H (^) (1 ■ 32+Tr) 1.01 0.78 

Relativistic three quark model §1 (z^l) ^'^^ 135 °- 68 

IMF model £ exp ( -(0.7) 2 Sfcf ) • ° ,7V r lF ^ — 1.44 0.66 

Skyrme model in the large N c limit H 0.99 exp (-1.3(w - 1)) 1.30 0.63 

MIT bag model H [J*-)^ 5 *^ 2 - 35 °- 45 



( p 2 f ) 1 = 0.66 ± 0.19 from B -> Drp 

(p2 f ) 2 = 0.71 ± 0.11 from £ -> D*l~V (5.5) 

with the world average values for V cb being \V cb \ = 0.0394 ± 0.0050 and \V cb \ = 0.0387 ± 0.0031, 
respectively. The weighted average of the two mesonic slope parameters are then p 2 M = 0.70 ±0.10. 
This translates into an upper bound for the baryonic slope parameter p 2 B according to the spectator 
quark model bound fl5.4|) . One has 

(p|)max = 0.89 ±0.19. (5.6) 

We mention that very likely the error on this bound will be considerable reduced in the near future 
with the new data expected from the bottom factories at SLAC and KEK. Combining both limits, 
( |5.3|) and (|5.6p , we obtain a prediction for the allowed values of the baryon slope parameter given 
by 

0.36 < p 2 B < 0.89 ±0.19. (5.7) 

According to these upper and lower bounds the first model (as remarked on before) and the last 
four models in Table I have to be excluded since they possess form factors which are too flat or too 
steep, respectively. The two QCD sum rule calculations [13,21] as well as the simple quark model 



evaluation [|^] feature slope parameters that satisfy the bounds fl5.7|) . We consider the two QCD 
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sum rule calculations to be the most reliable of the three model calculations since they are the 
least model dependent. Our final prediction for the range of values of the exclusive/inclusive ratio 
will be based on the two slope parameter values p 2 B = 0.85 and 0.65, resulting from the analysis 
of the diagonal and nondiagonal [^1] sum rules in fl21| , respectively. This range also includes the 
sum rule result of [13]. In determining our prediction for the range of ^ we shall also allow for 
the smaller inclusive rate calculated by the method of |[20]| . Thus our final prediction for the range 
of the exclusive/inclusive ratio in semileptonic — > A c decays is Re = 0.81 -4- 0.92. This range 
is consistent with the range of values from the semiquantitave analysis performed in Sec. II. Our 
conclusion is that the exclusive/inclusive ratio of semileptonic A^-decays is considerable higher 
than in the corresponding bottom meson case. 



VI. MISSING FINAL STATES 

Besides the quasi-elastic A5 — > A c contribution discussed before there are also A** resonant 
states and multiparticle final states contributing to the fully inclusive semileptonic A b rate. Of 
course, if the quasi-elastic contribution dominates the total inclusive rate much more than by the 
66% in the heavy meson case, there would not be much room left for the resonant and multiparticle 
final states. In the main body of this paper we have collected together theoretical evidence that the 
latter situation is very likely the case. One could turn this statement around in the following sense: 
if one would have theoretical reasons to believe that resonant and multiparticle final states are 
suppressed in inclusive semileptonic A^ — > X c transitions then the quasi-elastic exclusive A^ — > A c 
contribution must dominate. As we shall see there are theoretical reasons to believe in such a 
suppression in as much as some of the transitions to orbitally excited A** charm baryon states 
involve spin-orbit coupling transitions which are believed to be suppressed. 

The purpose of this section is to classify those final states in semileptonic A b — > X c transitions 
that form the complement of the quasi-elastic A;, — > A c transition. We divide these into class A 
contributions A5 — > A c Xlu, where the charm quark of the decay ends up in a charm A c directly 
or indirectly, and class B contributions, where the charm quark goes into a meson or a charm- 
strangeness baryon A^ — > X c (non A c )li>. Accordingly we define the two ratios 

T(A fe -> A c Xlu) 

RA - T(A b -> XM ■ (6 ' 1} 

and 



T{A b ^X c {nonA c )lu) 

Rb ~ T(A b XJu) ■ (6 - 2) 

Together with the exclusive/inclusive ratio defined before 

r(A b -> a» 

Re - r(A^x»' (6 ' 3) 

the three ratios must add up to one, i.e. 

R E + R A + R B = 1. (6.4) 



Note that all three ratios are positive definite which makes the constraint ( |6.4| ) potentially quite 
powerful if Re is close to one as is indicated by our analysis of the quasielastic rate in the previous 
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sections. As concerns the sizes of Ra and Rb one cannot even hope to provide semiquantitative 
answers at present. It is nevertheless useful to enumerate the final states belonging to the class A 
and class B transitions which we shall do in the following. 

Class A final states 



Potentially prominent among the class A final states are the transitions into the seven excited 
P-wave A**-states. Taking the bottom meson case for comparison theoretical estimates show that 
the corresponding transitions into excited mesonic P-wave states make up approximately 10% 
of semileptonic B-decays ^7|. The A**-states eventually decay down to the A c ground state via 
(multiple) pion emission or, with a much smaller branching fraction, via photon emission. There 
are altogether seven such P-wave states which are grouped into the three HQS doublets {A c ki}, 
{A c fci}, {A c fc 2 }, and the singlet {A c ko}- We use the terminology of [@] such that the excited K- 
and /c-states are symmetric and antisymmetric under the exchange of the momenta of the light 
quarks. The five symmetric states {A c ko}, {A^i}, and {A c k02} are made from a heavy quark 
and a light spin-one diquark. A5 transitions into these five states involve spin-zero to spin-one 
light-side transitions which can be expected to be strongly suppressed since they involve spin-orbit 
interactions. In the spectator quark model, where one neglects spin-orbit interactions, transitions 
into these five states are forbidden PEj| . It would be interesting to experimentally confirm this 
suppression. One thus remains with the transitions into the HQS doublet {A cK i} whose spin 1/2 - 
and spin 3/2~ members are very likely the recently discovered A c (2593) and A c (2625) states [|IJ. A& 
branching ratios into these states are not yet available. There could also be transitions into higher 
orbital A**-states. These transitions are, however, expected to be suppressed because of angular 
momentum suppression factors. Besides, transitions into symmetric higher orbital A**-states would 
again be suppressed due to spin-orbit coupling suppression. The suppression of transitions into 
the symmetric orbitally excited A**-states could be the source of the possible depletion of class A 
final states. For example, using spin counting only 1/3 of the existing P-wave excitations can be 
reached in semileptonic Ab transitions if the spin-orbit coupling suppression is active. 



d d d u c 



u b 




d u u u c 



u b 




a) b) 

FIG. 3. Class A final states. 

Another source of class A final states is accessible due to the creation of one or more additional 
(dd)- or (ml)-quark pairs in the basic transition. The relevant transitions for (dd) creation are (see 
Fig.3a) 

A°^A c + (£ c + )+X° f + r+^ (6.5) 
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or, when exchanging the d «-» u lines originating from the A&, one has 



A 



Xtf + l~+Vi. 



(6.6) 



For (mm) creation shown in Fig. 3b one has 



A 



X M + l +Vi. 



(6.7) 



The exchange of the d, u lines originating from the A& brings one back to ( |6.5D . Xm stands for 
a charmless mesonic inclusive state. Excited charm baryon states such as A** and X** are not 
explicitly included in the listing ( |6?5| - |6~7] ), but are implied. The S°, E+, and S+ + appearing 
in ( |6.5[ - |6.7| ) cascade down to the A+ state via pion emission making these processes class A final 
states. 

Class B final states 

There are two sources for class B final states. First there is (ss) quark pair creation where the 
strange quark ends up in a charm-strangeness baryon (Fig.4a) which decays weakly into noncharm 
states and therefore does not contribute to the class A final states []. Second, the charm quark of 
the decay may end up in a charm meson accompanied by (mm)-, (dd), and (ss)-quark pair creation 
as shown in Figs.4b-d. Let us list a few examples of such transitions. From (ss) pair creation one 
has (Fig.4a) 



A+ -> 5+ + X° + l- 



(6. 



or, when exchanging the d <-> u lines, one has 



a; 



(6.9) 



Xm 3 now stands for a strangeness meson state. Then there are the transitions where the charm 
quark goes into a charm meson. These are 



d s s u c 



d u d d c 



d u u u c 




d u b 




d u 





a) 



b) c) 

FIG. 4. Class B final states. 



d) 



3 The weak decay H c — > A c + n, though interesting, occurs only at the per mill level [29| 
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A+^D + + X° B + r+v h 



(6.10) 



A +^ D o + x + + r+Vh 



(6.11) 



(6.12) 



Xb stands for a light baryon state and Xb s for a strangeness baryon state. Excitations of the 
charm meson and charm baryon states are again implied. We do not discuss (cc) pair creation. 
The corresponding final states are barely accessible in semileptonic A& decays for kinematical 
reasons and will have a spectacular signature anyhow. 



We have brought together various pieces of theoretical evidence that the exclusive/inclusive 
ratio Re in semileptonic A b -decays is larger than in semileptonic B decays, where the exclu- 
sive/inclusive ratio amounts to 66%. We predict that the exclusive quasielastic semileptonic A 6 - 
decays make up between 81% and 92% of the total inclusive semileptonic A b rate. At present 
there is no experimental information on either the exclusive or the inclusive branching ratio in 
semileptonic Ab decays. The problem is that present and planned experiments do not have access 
to reliable A b tags which are necessary for a measurement of their branching fractions. Ideally 
one would run a e + e~-machine right above A;,A fe threshold which would solve the tagging problem. 
However, such experiments are not planned in the foreseeable future. The above assertion about 
the dominance of the quasi-elastic mode in semileptonic A b decays may take a long time to verify 
experimentally. It may nevertheless be used as a working hypothesis in the experimental analysis 
of semileptonic A& decays in particular if further theoretical progress in the theoretical description 
of semileptonic A& decays confirms the estimates made in this paper. 



This investigation was prompted by two questions of our experimental colleagues P. Roudeau 
and G. Sciolla who asked us about theoretical expectations for the size of the the A c + X (class 
A) and (nonA c ) + X c (class B) contributions in inclusive semileptonic A& decays. We did our best 
to try and provide answers to these questions in our paper, at least partially. We would like to 
thank N. Uraltsev for fruitful discussions. We also thank O. Yakovlev who participated in the early 
stages of this calculation. The work of B.M. was partially supported by the Ministry of Science 
and Technology of the Republic of Croatia under the contract Nr. 00980102. 



VII. SUMMARY AND CONCLUSION 
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